Tailoring symmetry groups using external alternate fields. 
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Macroscopic systems with continuous symmetries subjected to oscillatory fields have phases and 
transitions that are qualitatively different from their equilibrium ones. Depending on the amplitude 
and frequency of the fields applied, Heisenberg ferromagnets can become XY or Ising-like — or, 
conversely, anisotropies can be compensated — thus changing the nature of the ordered phase and 
the topology of defects. The phenomenon can be viewed as a dynamic form of order by disorder. 

PACS numbers: Q5.50.+g, 75.40 G. 
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The effect of alternate fields on interacting many-body 
systems has been a subject of long-standing interest. In 
particular, magnetic hysteresis has been extensively stud- 
ied for clean or disordered Ising-like systems [1-5] . 

A much less studied problem is the effect of a periodic 
drive on systems having a continuous symmetry either in 
a pure or a disordered [7-9,6] case; even though, as we 
shall see, the existence of soft modes make the systems re- 
spond in striking and interesting ways. Years ago, Rao et 
al. [7,8] and Dhar and Thomas [9] studied a driven clean 
A^-component Heisenberg ferromagnet, solvable analyti- 
cally in the large- A^ limit for any dimension d. Due to the 
fact that at N = oo almost any configuration is perpen- 
dicular to the field, one can not study thus, phenomena 
that depend on a delicate balance between a tendency of 
the system to order longitudinally or transversely to the 
field. 

In this paper we explore the full problem for d > 2 and 
A^ > 2. We firstly treat the XY {N = 2) model and show 
that, surprisingly enough, there are in fact three kinds 
of ferromagnetic phases, with the magnetization longi- 
tudinal, transverse and canted with respect to the field, 
respectively. We show this directly for the mean-field 
{d ~ oo) case, and extend the result for all dimensions 
d > 2 by use of low frequency and low temperature ex- 
pansions. The same analytic methods show that for the 
Heisenberg case (A^ > 3) there is always transverse order 
[9], again in all d > 2. In addition we show that one can 
tailor in this case the symmetry group from 0(3) to 0{2) 
or to Z2 by applying one or more a.c. fields: a possibil- 
ity we argue is quite general for systems with continuous 
symmetries. 

Our motivation is that an understanding of this sys- 
tem is a base to explore the connections with wider class 
of forced systems with continuous symmetries, including 
liquid crystals [10], lamellar polymers [11], crystalline de- 
fects, ferromagnetic conductors [12], etc. 

Consider the 0{n) ferromagnet 
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where 2d is the number of neighbors, itself depending 
on lattice dimensionality and topology. The S are A^- 
dimensional vectors either of fixed norm (hard spins ), or 
whose length may fluctuate according to a soft-spin term. 
The system is coupled to an ac field: 



Eh = — cos{ujt) h ■ 



(2) 



and we shall consider the dynamics in the strongly dissi- 
pative limit: 



ST = - 



dE 



(3) 



where r^f is a white gaussian noise of variance 2T. The 
energy E is the sum of (1), (2) and a soft-spin term fixing 
the spin length or, alternatively, a Lagrange multiplier for 
hard spins. In (3), we are neglecting precessional effects 
(which bring in a host of interesting phenomena [13]). 

Selection: order by disorder. 

Consider the 'hard' spin model at zero temperature, 
subjected to an a.c. field. We shall use the example of 
the XY model, but the argument is valid for any number 
of components and in any number of dimensions. Using 
polar coordinates, we have for the angle 0^: 



9j) ~ h cos((jji) cos ( 



(4) 



This always admits a solution in which all spins move in 
phase 9i — 9, Vi (inhomogeneous solutions decay into 
this one) which makes interaction terms zero and we have 
for any geometry: 



9{t) = 2 tan" 



(5) 



By considering all possible values of the integration con- 
stant k, we conclude that at zero temperature solutions 
are possible in which the total magnetization vector os- 
cillates around any possible angle, and without hystere- 
sis. Unlike the case without field, these solutions are no 
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longer related by a continuous symmetry, and only the 
discrete symmetries 9 —9 and 9 ^ n — 9 remain. (In 
the case of A'' > 2, a symmetry subgroup survives: for 
example, in A'^ = 3 the spins evolve with 9 as in (5) and 
with constant (p; solutions with the same k but different 
if are related by an 0(2) symmetry.) 

In situations like this one, when one has at T = solu- 
tions that are unrelated by a symmetry, thermal or quan- 
tal fluctuations will generically select a subset of them 
(unless of course they destroy the order altogether) . This 
phenomenon has been named order by disorder [14], as 
it is the fluctuations that are responsible for a reduction 
in the number of solutions. In practice, if we prepare 
the system at small T in one of the possible T = solu- 
tions, it will drift to the selected average angle: there is 
a secular perturbation acting on timescales much slower 
than the vibration frequency. Temperature also brings in 
hysteresis and hence dissipation. 

Paramagnetic, transverse, longitudinal and 
canted solutions. 

In an 0{N) ferromagnet, one can have four types 
of phases, which can be best distinguished by consid- 
ering the magnetization in the direction parallel and 
perpendicular to the field Mh{t) and M±{t), the angle 
9{t) = tan^ ^{M±{t)/Mh{t)), and the corresponding aver- 
aged quantities: M« = ^ f dtM°'{t) &nd9=\§dt 9{t) 

• Paramagnetic: The magnetization follows the fleld 
with a delay (hysteresis): M±{t) = and Mh = 0. 

• Longitudinal {9{t) = or 9{t) = n): the 
magnetization points in the direction of the field: 
Mj_(i) = 0, but Mfc ^ 0. 

• Transverse {6 = Tr/2 or 9 = —Tr/2): the magnetiza- 
tion has a non-zero component M±{t) ^ orthogo- 
nal to the field, the component parallel to the field 
has zero time-average Mh = 0. 

• Canted {0 < 9 < 7r/2 or Tr/2 < 9 < it): The mag- 
netization evolves around an oblique angle with the 
field's direction: Ml ^0 undM^^O. 

The longitudinal solution is doubly degenerate. In the 
XY case, the canted solutions have degeneracy four and 
the transverse ones two. Both of them are continuously 
degenerate for > 2 (one solution per plane determined 
by iMh,M^)). 

Dynamical phase transitions to a magnetized phase arc 
well attested in the case of an Ising system [15-17]. For 
continuous 0{N = oo) systems, Rao et al ( [7,8]) found a 
dynamical phase transition on increasing frequency from 
a paramagnetic state with M = to a ferromagnetic 
regime. Subsequently, Dhar and Thomas [9] pointed out 
that the order is in fact always transverse, and Mh = 
[18]. They also studied the case > 2 in finite dimen- 
sions within the unmagnetized phase [19,20]. 



Dynamic mean-field approximation. 

The dynamic mean-field approximation has the advan- 
tage of being completely solvable, and one can easily get 
a complete phase diagram. The equations consist of a 
single-spin equation and a self-consistency condition: 

S-" = (M" + h°'){t) - X{t)S" + rf 

M"(<) = (5"(^))s^„g^esp^„ (6) 

where A is a Lagrange multiplier imposing the spin 
length. One can solve the system (6) by considering the 
evolution of the set of expectation values. For example, 
for the XY model, using the Fokker- Plank equation as- 
sociated with (6) we write an exact infinite system of 
equations for a;„ = (cos(n^)), ?/„ = (sin(n^)), n = 1, ...: 

Vn = -n^Tyn + ^nxi{yn-i - y„+i) + m(x„_i + x„+i) 

in = -n^Txn + ^nxi{xn-i - Xn+i) - m(2/„_i + y^+i) 

yi = - T)yi - ^x^y2 -\{Vi + Kt))x2 + \h{t) 

a;i = (^ - T)xi - ^xrX2 - ^(yi + h{t))y2 (7) 

with m{t) ~ n{yi + h{t))/2. For the Heisenberg case one 
has to study the dynamics of the expectation values of 
the spherical harmonics Yim- 

We have numerically solved the system (7), keeping as 
many modes as necessary, for various values of h,w,T. 
The main results are summarized in the figures. 
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FIG. 1. ^ as a function of uj for the XY mean field model, 
h = 1, T = 0.2. The transitions are second order. 

FiG. 1. shows the value of ^? vs a; at temperature 
T = 0.2 and h = 1. We see second order transitions, 
as frequency decreases, from longitudinal (L) to canted 
(C), and from canted to transverse (T). FIG. 2. gives the 
general phase diagram of the system at /i = 1 (the dashed 
line at T=0.5 represents the para-ferro transition at zero 
field). In the limit of high frequencies lu oo, h/uj 0, 
there is a dynamical critical temperature ~ .42 from 
a longitudinal to a transverse phase. 
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FIG. 2. T-uj Phase diagram for the XY mean field model, 
h = l. 



In the Heisenberg case, we have only found paramag- 
netic or transverse behavior. Below we confirm this by 
showing that, in all dimensions d > 3, no longitudinal 
order appears in the first order in low-temperature ex- 
pansion (the most favourable case, as seen in Fig. 2), 
thus supporting the conjecture of [9] that the Heisenberg 
model, and a foriori all the N > 2 cases, are in this sense 
qualitatively different from the XY model. However, as 
we shall see below, one can still select a single direction in 
a Heisenberg case, but with the application of two fields. 

Let us stress here that at very low frequencies mean- 
field is potentially very misleading, as the physics may 
be dominated by nucleation effects which this approxi- 
mation neglects. 

Weak fields cind low frequencies: a general 
mechanism for transverse selection. 

In the limit of low frequencies and weak fields (w ^ 0, 
h/w constant), a general simple argument can be given 
to show that there is transverse selection for any > 3 
and d > 3: Consider a field applied along the z axis, 
and split it in components aligned hM{t) and normal 
h±{t) to the instantaneous magnetization. Because by 
assumption the field is weak and varies slowly, the ef- 
fect of the aligned component of the field is to change 
the magnetization norm M linearly and adiabatically: 
M{t) = Mo + Chuit), where Mq is the unperturbed 
norm and C a positive susceptibility. On the other hand, 
the effect of the perpendicular field is to modify the an- 
gle between the magnetization vector and the equato- 
rial plane as: J{M)Cl = h±{t). The mobility coefficient 
J^^{M) depends only upon the norm M due to symme- 
try. Putting this together we get: 

J{M)n = h{t) cos n ; M = Mo + Ch{t)smQ (8) 

and using the weak field approximation: 

//(O cosi). 



It is easy to check that this equation yields transverse 
selection provided J'(Mo) > 0, i.e. when the mobil- 
ity decreases with the magnetization. Let us note that 
an inhomogeneous nucleation mechanism would need a 
critical droplet size of h^^^^, and a free-energy barrier 
~ which in the present regime requires much 

longer times compared to uj~^. 

We can now explain in words the selection as a ratchet 
mechanism: suppose the magnetization starts at an an- 
gle in the first quadrant. During the semicycle when the 
magnetic field points upwards it simultaneously rotates 
the magnetization towards the z axis and it stretches 
the norm. In the negative cycle, the rotation is towards 
the xy plane and the norm is shortened. The slight 
changes in the norm make the positive cycle less efficient 
(if J' (Mo) > 0) than the negative cycles, and hence the 
vector has a net drift towards the transverse direction on 
each cycle. The same mechanism applies for soft spins in 
all dimensions at zero temperature. 

We have checked that indeed J' {Mo) is positive in all 
cases. Near the critical point J(M) oc M, implying that 
the selection mechanism dissappears at the paramagnetic 
transition. 

Low T expansion: longitudinal selection 

These models arc obviously exactly solvable at T = 0, 
and an expansion in powers of T is feasible in any dimen- 
sion. In the XY case, one decomposes 9i{t) = {Oi{t)) +9i 
(averages over the thermal noise), and one writes an 
evolution equation for the expectation values (9i){t), 
{()iej){t), {ei0j0k){t), .... To leading order in T only two- 
point correlations are necessary. In the most general case, 
we have: 



^ = h{t){i-l{el))cos{e.) 



9,))-/i(Mfe)sin(0„) 
(10) 



n 



J{Mo) + CJ'{Mo)h{t) sin f2 



(9) 



where Aai oc g if a, i are neighbors in d-dimensional space 
and zero otherwise. This system is linear in the corre- 
lations. One can assume translational invariance with 
one-point functions independent of the site and two-point 
functions depending only on the distance between sites. 
The system then becomes exactly solvable in Fourier ba- 
sis (details in [21]). For hard spin systems, it always 
yields longitudinal selection for any d > 3. Thus, we 
have shown the existence of a longitudinal and a trans- 
verse phase for the XY- model in finite dimension d> 3. 
The existence of the intermediate canted phase we have 
obtained analytically within mean-field (large d) has been 
checked with simulations. 

The generalization to the Heisenberg {N = 3) case of 
the low-temperature expansion is immediate, one has to 
consider two-point correlations of both angles 9i and 0, 
[21]. One can show that there is no selection to order T 
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in any dimension d > 3. Since the lower temperatures 
are the most favourable for a longitudinal order, this is 
evidence that there is no longitudinal phase at all. 

High frequencies 

Another instructive method is the high-frequency ex- 
pansion in which the external field can be treated pertur- 
batively. We start from a system in equilibrium under an 
infinitesimal field in a direction ^. In terms of the linear 
X^p and higher susceptibilities, the variation of the mag- 
netization vector Ma in presence of the a.c. field ha{t) 
is given to second order as: 

j dt'xip{t,t')hp{t') + j dt'dt''xifj^{t,t\f'')h8{t')h,{t'') 

The linear term does not contribute to the drift if 
f ha{t)dt = 0. The quadratic term does and, in the large 
CO limit, the contribution reads: 



for these out of equilibrium problems, since there are few 
qualitative ideas to guide us before having an actual solu- 
tion, as we have for equilibrium thermodynamic systems. 
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lim 



dt' 



(11) 

where hphj = § hij{t)h^{t)dt ^ 0. Specialising to the 
mean-field case, already the first order in T and in 
gives for the XY case a longitudinal selection correspond- 
ing to an effective potential ~ Th^ui~^ sin(20). 

Now, if the field is a sum of two components acting at 
right angles with a ratio of frequencies larger than two, 
interferences appear to the third order and their contri- 
bution is additive up to w"^. We have checked in the 
hard spin case that two orthogonal fields at two high 
enough frequencies oj and Sw, make the Heisenberg sys- 
tem Ising-like: the magnetization is transverse to one 
field and longitudinal to the other. The same effect can 
be obtained with a field with constant modulus that ro- 
tates on a plane. 

Perspectives 

One can forsec several applications of these ideas in 
different fields. Changing the symmetry group of a sys- 
tem leads to a change in the topology of the defects: 
thus an alternate field may perhaps be used as a tool to 
study in detail the role these play in the transitions [12], 
in the phase-ordering kinetics, in the anomalous Hall ef- 
fect (since it involves the interaction of electrons with 
the topological defects) [12], etc. Similarly, the chiral- 
ity scenario [22] for the spin-glass transition (which as- 
sumes that spin glasses are essentially isotropic) could 
be put to test experimentally by monitoring the effect of 
an increased and a decreased effective anisotropy. Sub- 
jecting ordered liquid crystals to a.c. fields might be a 
way to directly observe what becomes of the order and 
the textures, a strategy for which there are interesting 
precedents [23]. 

Such richness of amusing and perhaps useful phenom- 
ena is in contrast with the extreme poverty of concepts 
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